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ABSTRACT 

This paper presents results of a survey search for bright compact radio sources at 22 GHz with the 
VERA radio-interferometer. Each source from a hst of 2494 objects was observed in one scan for 
2 minutes. The purpose of this survey was to find compact extragalactic sources bright enough at 
22 GHz to be useful as phase calibrators. Observed sources were either a) within 6 degrees of the 
Galactic plane, or b) within 11 degrees from the Galactic center; or c) within 2 degrees from known 
water masers. Among the observed sources, 549 were detected, including 180 extragalactic objects 
which were not previously observed with the very long baseline interferometry technique. Estimates 
of the correlated flux densities of the detected sources are presented. 
Subject headings: radiosources — VLBI — catalogues — surveys 



1. INTRODUCTION 

VLBI Exploration of Radio Astrometry (hereafter 
VERA) is a Japanese VLBI array dedicated to phase- 
referencing astrometry to explore the Milky Way Galaxy 
(|Honma et al. 2000a; Kobavashi ct al. 2004). VERA ob- 
serves Galactic H2O and SiO maser sources and deter- 
mines their parallax and proper motions at the 10 mi- 
croarcsec accuracy level. Based on the high-precision 
astrometry of ~1000 maser sources, VERA will study 
the three-dimensional structure and dynamics of the 
Galaxy's disk and bulge, revealing the true shape of the 
bulge and spiral arms, its precise rotation curve and the 
distribution of dark matter. 

The VERA array consists of four 20-m diameter radio 
telescopes spread over Japan, covering baselines ranging 
from 1000 km to 2300 km. Unlike other VLBI stations, 
VER A telescopes are equip ped with a dual-beam sys- 
tem ijKawasfuchi et al1l200(t >. with which one can simul- 
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taneously observe two sources separated from each other 
by up to a hundred times the telescope's primary beam 
size. With such a system, VERA can observe a Galac- 
tic maser source and an extra-galactic calibrator source 
at the same time and perform relative astrometry of the 
Galactic maser source with respect to the extra-galactic 
calibrator, which serves as a position reference. Simulta- 
neous observations of two sources, rather than switching 
observations with single-beam telescopes, are expected to 
be more effective in canceling out the tropospheric fluctu- 
ations, which are the main source of position error in dif- 
ferential VLBI astrometry at 22 GHz or higher frequency. 
In fact, preliminary results from VERA's dual-beam sys- 
tem observations have already shown its high capability 
of phase-refere ncing to reach down to 10 microarcsecond 
level accuracy ()Honma et al.ll200'3l) , and now monitoring 
of several bright maser sources is on-going to produce 
initial results from the VERA. 

In order to observe a pair of two adjacent sources with 
VERA, the separation angle of the pair must be less than 
2.°2. This is not only a mechanical limit of VERA's 
dual-beam system, but also a practical limit for effec- 
tive phase referencing to achieve 10 microarcsec accu- 
racy, i.e., the reference source needs to be close enough 
to the target maser to remove the tropospheric fluctua- 
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Fig. 1. — VERA network. 



tions sufficiently. Therefore, the existence of calibrators 
close to target sources is crucial to successful VERA ob- 
servations. If we assume a uniform density of calibrators 
in the sky plane, ^ 2700 calibrators are required in the 
whole sky to find a calibrator within 2.° 2 of any sky po- 
sition (a circle with radius of 2.°2 covers 0.0046 steradi- 
ans in the sky; thus 47r/0.0046 w 2700 are necessary). 
Practically though, more sources are needed because the 
distributions over the sky is not uniform. 

To date, there have been many efforts to find good 
ca librators observab l e with VLBI. The pioneering work 
of iCohen fc Shaffer! l)1971[) provided the first catalogue 
of source positions determined with VLBI, which con- 
tained 35 strong sources. In 1998 the International Ce- 
lestial Reference Frame (ICRF) catalogue of 667 sources 
produced by analyzing VLBI observations made in the 
framewor k of sp ace geodesy programs was published by 
iMa et all l|1998() . At present, the most massive survey 
of compact calibrators is the VLBA Cali brator Survey, 
hereafter VCS, performed in 1994-2005 fBeaslev et al.' 
^002; Fomalont et al. 2003; Petrov et al. 2005, 2006), 
and Kovalev et al. (submitted to AJ, 2006^). After the 
fifth VCS campaign, the NASA catalogue 2005f_astro^ of 
all sources detected in the absolute astrometry mode at 
2.2 GHz and/or 8.6 GHz frequency bands (S and X) con- 
tains 3481 sources. Among them, 85% are acceptable as 
calibrators for the VLBA. Nevertheless, in the sky area 
at 6 > —40°, which is observable with the VERA, the 
probability of finding an X/S calibrator with correlated 
flux density over 60 mJy within 2.°2 of any position is 
~70% (Petrov et al. 2006), indicating that more calibra- 
tors are needed in order to observe all maser sources with 
VERA. Moreover, maser sources tend to be confined to 

^ Available at http : //arxiv ■ org/ast ro-ph/0607524] 
^ Available at |http : //vlbi ■ gsf c ■ nasa . gov/astro J 



the Galactic plane, because they are basically located in 
the Galaxy's disk, but the density of known calibrators 
near the Galactic plane is even lower than that in other 
sky regions. This is due to the strategies used to find 
compact calibrators: the list of candidates was based on 
low resolution sky surveys that severely suffered from 
strong radio foreground in the Galactic p lane. Besides, 
some past surveys, JVAS llWrobel et alJll998|) . VCSl, 
VCS3, VCS4, specifically avoided the Galactic plane. To 
overcome this problem, a calibr ator survey in th e Galac- 
tic plane has been done bv .Honma et alJ l)2000bf) finding 
more than 50 new calibrators at |&| < 5°, but still many 
more calibrators in this region are needed. 

For these reasons, we have been performing the Galac- 
tic Plane Survey of compact radio sources, hereafter 
GaPS, based on the VERA and VLBA observations. The 
objective of this campaign is to provide a list of calibra- 
tors for dual-beam VERA observations of H2O masers 
at 22 GHz. Our strategy for GaPS is two-fold: first, we 
observe a wide list of candidate sources with the VERA, 
and then observe a narrow list of detected sources with 
the VLBA in order to determine their positions at the 
1-10 nrad level of accuracy and to obtain high quality 
maps. In this paper, we present results from the first part 
of the GaPS project, the fringe search observations with 
the VERA, reporting the detection of 180 new compact 
sources. In section ^ we describe our approach to the 
selection of candidates, scheduling and observations. In 
section |21we describe analysis of the observations. Re- 
sults are discussed in section SQlE^J^d concluding remarks 
are given in section [JS] 

2. OBSERVATIONS 
2.1. Source selection 

We searched for compact extragalactic sources with 
correlated flux density > 100 mJy at spacings ~ 10^ 
wavelengths at 22 GHz in the following zones: a) |6| < 6° 
and S > —40°; b) within an 11° radius of the Galac- 
tic center; and c) within 2.°2 of objects listed in the 
2nd update of Arcetri Catalog of H2O maser sources 
iValdettaro et al] |200lD and the stellar masers catalogue 
of lBenson et alJ l|199nD . 

First, we combed through the wide list of 3481 known 
compact sources observed under absolute astrometry 
programs at S/X bands. Among them, 510 fall in our 
zones of interest. Also among the 3481 known compact 
sources, 252 objects have been observed u nder the K/Q 
band survey in 2002-2005 l)Fev et al.ll2005t I Jacobs et alJ 
(2005), and can be considered as confirmed K-band cali- 
brators. But only 16 of these fall in our zones of interest. 
The first list of targets was formed from the 510 known 
X-band calibrators, except 16 objects observed with the 
K/Q band survey. The purpose of including known X/S 
calibrators in the list of candidate sources was to check 
whether they are bright enough at K band to be detected. 
It was expected that not all of these will be detected with 
the VERA, since although they have correlated fiux den- 
sities above 60 mJy at X-band, their correlated flux den- 
sities at K-band may be below the VERA detection limit, 
and, therefore, not suitable as calibrators for VERA. 

Second, we searche d the CATS database 
ijVerkhodanov et al.l 11997(1 containing almost all ra- 
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dio catalogs'^ to find flux density measurements at radio 
frequencies 1.4 MHz and higher. We selected the sources 
with at least two measurements of the flux density 
and found their spectral index by fitting a straight line 
to the log(F)/log(f) dependence, where F is the flux 
density and / is the frequency. The flux density was 
extrapolated to 22 GHz. Those sources in the zone 
of our interest, which had extrapolated flux density 
> 100 mJy at 22 GHz, and which were not observed in 
the absolute astrometry mode, were put in the second 
list of targets, in total 1901 new sources. 

2.2. Scheduling 

The list of target sources was split into two group: 
the first group, of 1386 objects, within 6° of the Galac- 
tic plane or within 11° from the Galactic center; and a 
second group, of 1032 sources, within 2.°2 from known 
galactic maser sources. The first group had priority in 
scheduling. Among target sources, 94% were observed in 
one scan of 120 s, and 6% were observed in two scans. 
Observing schedules were prepared using the following 
procedure. A source in the target list, which was not ob- 
served in previous experiments and which had the highest 
elevation at the station of VERZMZSW, was selected as 
the first objects in the schedule. In order to select the 
next source, the pool of target sources was examined, and 
for each source, which was at least 10° above the horizon 
and which was not previously observed, a score was com- 
puted: S = 100/Ts;e^ +30* (20 -S) + 0.2 * Emin where 
Tsiew is the slewing time in seconds, 6 is the source decli- 
nation, in degrees, and Emin is the minimum elevation of 
the source among the network stations, in degrees. The 
source from the first group with the maximum score was 
put in the schedule. If no source from the first group was 
visible, the source with the highest score from the second 
group was put in the schedule. Then the procedure was 
repeated. The slewing rate of VERA stations is 2° 
and acceleration is 2° s~^. The antenna control hardware 
requires adding a margin of 10 seconds to slewing time, 
and the slewing time should be rounded to the next 4 s. 
Measurement of system temperature takes an additional 
30 seconds. System temperature was measured either ev- 
ery 5 scans or when the antennas moved more than 10° 
from the previous measurement, whichever occurs first. 
In addition to targets, observations of amplitude calibra- 
tor sources were scheduled every hour. At the beginning 
and at the end of each session a very strong source, a 
fringe finder, was observed. The amplitude calibrator 
sources were taken from the list of 252 sources observed 
with the K/Q band survey. The average density of ob- 
servations was 22 scans per hour. 

2.3. Observations 

GaPS survey observations were scheduled mainly 
within gaps between regular VERA observations (ta- 
ble ^ . Some observations were performed only at three 
antennas while one of the stations was not available. 

The left circular polarization in band 21.97-22.47 GHz 
band was received, sampled with 2 bi t quantization, and 
filtered using the VERA digital filter ijiguchi et al.ll2005j) 
before being recorded onto magnetic tapes. The digital 
filter split the data within the 500 MHz band into 16 
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TABLE 1 
Observing sessions 



(1) 


(2) 






(3) 


(4) 


(5) 


r05284b 


2005 Oct. 


11 


08:20 


7.0 


3 


1.0 


r05285b 


2005 Oct. 


12 


08:20 


7.0 


4 


1.0 


r05332b,c 


2005 Nov. 


28 


12:20 


10.8 


4 


1.0 


rUoooob 


2000 Dec. 


01 


11:00 


12.5 


4 


1.0 


r05343b 


2005 Dec. 


09 


15:00 


8.5 


4 


0.1 


r05348a 


2005 Dec. 


14 


09:00 


14.0 


3 


1.0 


r05349a 


2005 Dec. 


15 


10:00 


20.0 


3 


1.0 


r05350a 


2005 Dec. 


16 


07:00 


7.5 


3 


1.0 


r05350b 


2005 Dec. 


16 


14:55 


8.5 


3 


1.0 


r05353d 


2005 Dec. 


19 


04:00 


5.5 


3 


1.0 


r05359d 


2005 Dec. 


25 


02:00 


7.0 


4 


0.1 


r06052b 


2006 Feb. 


21 


18:00 


13.0 


4 


0.1 


r06054a 


2006 Feb. 


23 


19:00 


11.0 


4 


0.1 


r06059b 


2006 Feb. 


28 


18:00 


16.0 


4 


0.1 


Note. — 


(1) session name; (2) start date; 


(3) 



duration in hours; (4) # stations; (5) Accumulation 
period length in seconds. 



frequency channels of 16 MHz width each, equally spaced 
with spacing 32 MHz. The data were recorded with the 
SONY DIR2000 recorder at 1 Gbps rate. 

3. DATA ANALYSIS 

Correlation was carried out on the Mitaka FX corre- 
lator (Chikada et al. 1991!) with a spectral resolution of 
250 kHz. Thus, the amplitude and phase of the cross 
correlation function at 64 spectral channels within each 
of the 16 frequency channels was determined. The corre- 
lator writes the output in the native CODA format. The 
output was re-formated to the FITS-IDI format (C. Flat- 
ter, (1998) unpublished^). One part of the observing ses- 
sions was correlated with accumulation periods of 0.998 s, 
and another part was correlated with accumulation pe- 

rinrlQ nf 1 q 




Stl: VERAIRIK St2: VERAOGSW Sou: 0B35-339 Exp: r0533Kb-nol Step: 2.424240-10 s 4.347830-14 a/t 

Obs: 602 Scan: 332-1654 Start: 2005.11.36-16:54:40.563 DuraUon: 119.931 sec 



Ampl: 0.0000690 SNR: 11.49 Processed on: 2006.07.30-22:46:38 

Gr_del; -7.507202D-09 sec Pll_rat; -4.112511D-13 s/s PIMA v 0.78 3005.07.11 

Fig. 2. — Fringe amplitude versus delay and delay rate for a 
strong source. 



The data were analyzed with program PIMA. The 
first step of analysis was fringe fitting. The group de- 
lay and delay rate that maximize the fringe amplitude 

^ Available on the Web at 

Ihttp : //www . c V . nr ao . edu/f its/documents/draf ts/idi-f ormat .ps 
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Plotting window: 

Delay window [ -9.11059D-08, -1.1 1059D-0B] a 
Rate window: [ -2.26710D-ia. 1.73290D-12] s/a 
Stl: VERAMZSW 312: VERAIRIK Sou: 2034+497 Exp: r05335b-nol Step: 3.225810-10 s 4.347830-14 a/a 

Obs: 171 Scan: 335-1218 Start: 2005.12.01-13:18:33.123 DuraUon: 118.931 sec 

Ampl: 0.0000228 SNE: 5.62 Procesaed on: 2006.08.19-19:11:21 

Gr_del; -5.110586D-08 sec Ph_rat; -2.670961D-13 s/s PIIU v 0.313008.03.14 

Fig. 3. — Fringe amplitude versus delay and delay rate for a 
marginally detected source. 



were sought for each baseline independently. The spec- 
trum of the cross-correlation function was presented as 
a two-dimensional array, with the first dimension run- 
ning over frequency channels, and the second dimension 
running over accumulation periods. The size of the ma- 
trix was 960 X 120 for scans with accumulation periods 
of 1 s and 1920 x 1200 for data with accumulation peri- 
ods of 0.1 s. The grid was extended to 4096x512 and 
8192 x4096 for the data with 1 s and 0.1 s accumula- 
tion periods, respectively, with an oversampling factor 
of four by padding extended elements with zeros. The 
amplitudes of the cross-correlation function was normal- 
ized by dividing it by the square root of the product of 
the autocorrelation of the recorded data at each station 
of the baseline. The two-dimensional fast Fourier trans- 
form was performed. The first dimension of the Fourier 
transform runs over group delays from through 4 mks 
with steps of 0.49 ns, and the second dimension runs 
over delay rates from through 4.6 • 10^^^ with steps of 
0.9 • 10~^^ for the data with 1 s accumulation periods, 
and from through 4.6 • 10"^° with the step of 1.1 • lO^^^ 
for the data with the accumulation periods of 0.1 s. The 
coarse search for the maximum of the correlation ampli- 
tude was performed over results of the Fourier transform. 
In the worst case, when the maximum of the correlation 
amplitude falls just between the nodes of the grid, the 
amplitude is reduced by sinc^(7r/(2w)), where v is the 
oversampling factor. When w = 4, in the worst case the 

amplitude is reduced by ^^(^^^) — 0.95. The oversam- 
pling is important for detecting weak sources. The fine 
search of the correlation maximum is performed by the 
consecutive direct Fourier transform at a grid 3x3 around 
the maximum with reducing spacings of the grid by half 
at each step of an iteration. Three iterations of the fine 
search were performed. Examples of fringe amplitude 
plots versus trial delay and delay rate for a weak and a 
strong source are shown in figures El and |21 

Observations of strong sources with signal to noise ra- 
tios exceeding 120 were used for evaluation of a baseline 
dependent bandpass calibration. The phase of the band- 
pass calibration was determined as a residual phase with 
the opposite sign at each spectral channel after subtract- 



ing the multi-band fringe phase. A linear model was 
fitted into individual phases at spectral channels within 
each frequency channel. The amplitude of the band-pass 
was determined as the ratio of the fringe amplitude at an 
individual spectral channel to the maximum amplitude 
within the band. The estimate of the maximum ampli- 
tude across the band was found as the maximum value 
of the amplitude across 1024 spectral channels minus the 
doubled standard deviation of the amplitude at an indi- 
vidual spectral channel. It was found that the bandpass 
is stable between experiments, so the bandpasses over 
all experiments were stacked and averaged. An exam- 
ple of the bandpass calibration is shown in figures 01 and 
El The complex bandpass calibration was applied to all 
observations, i.e. the uv data were multiplied with the 
complex bandpass before the fringe search. This tech- 
nique increased the fringe amplitude by 5-15%. 



0.2 




Frequency MHz 

Fig. 4. — The phase of the bandpass calibration at baseline 
VERAMZSW / VERAOGS W. 




Frequency MHz 

Fig. 5. — The amplitude of the bandpass calibration at baseline 

VERAMZSW / VERAOGS W. 



3.1. Calibration 

System temperatures including atmospheric attenu- 
ation were measure d with the chopper-wheel method 
(lUlichfc HaasllT97l . A microwave absorber at ambient 
temperature was inserted just in front of the feed horn at 
the beginning of each scan, and the received total power 
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was nieasured with a power meter. Using the measured 
total power for the blank sky and the absorber, we can 
determine the temperature scale automatically corrected 
for the atmospheric attenuation. The system tempera- 
tures during the best weather conditions were 100-200 K. 
However, at typical weather conditions, T^y^ were 200- 
500 K depending on the elevation of observed sources, 
reaching 2000 K at low elevations and during adverse 
weather conditions. Figures El and [3 illustrate recorded 
system temperatures at two stations at good and adverse 
weather conditions. The uncertainty in the temperature 
scale was estimated to be 10%, mainly due to the as- 
sumption that the ambient temperature is the same as 
the air temperature 



2000 



1500 



1000 |5 



E-h" 500 - 




2x10 



4x10^ 



Time s 

Fig. 6. — System temperature corrected for the atmosphere opac- 
ity versus time from beginning of observing session at VERAISGK on 
2005.12.14. 



250 



200 



150 - 




100 - 



80 

Elevation deg 

Fig. 7. — System temperature corrected for the atmosphere opac- 
ity versus elevation in degrees at VERAMZSW on 2005.12.09. 

The aperture efficiency of each antenna was 45-52% 
in the 22 GHz band, which was measured by observing 
the continuum emission from Jupiter once per year. The 
accuracy of the aperture efficiency measurement was typ- 
ically 2-5%. The antenna gain does not show elevation 
dependency in the elevation range 10-90°. 

Amplitude calibration was performed in two steps. 
First, the amplitude conversion factor was computed 
for each observation based on measurements of system 

temperature and antenna efficiency as ^ = 



Then the multiplicative correction to the conversion fac- 
tor was computed for each observation of the amplitude 
calibrator at each baseline as SA = -p^ where Fkq is 
the correlated flux density from maps from the KQ sur- 
vey and Fobs is the observed correlated flux density. The 
correction SA was averaged over all amplitude calibra- 
tors within an individual observing session. The scat- 
ter of the individual estimates of the correction provided 
the measure of errors of amplitude calibration: 20-25%. 
Since for the vast majority of sources, only 3 or 6 val- 
ues of fringe amplitude and phase were determined, no 
attempt to make maps was made. 

3.2. Detection limit 

Strong sources have a noticeable peak in plots of fringe 
amplitude versus trial delay and delay rate, for instance 
figure[21 The presence of noise will force the fringe search 
process to find a peak in such a plot even in the absence of 
a signal from the source. For each observation we need to 
evaluate the probability of false detection. In the absence 
of signal the amplitude of the cross correlation function 
z has Rayleigh distribution (Papolius 1965.') : 



p{z) 



(1) 



where a is the standard deviation of the real and imag- 
inary part of the cross correlation function. The noise 
level of fringe amplitudes was evaluated for each obser- 
vation. At the grid of cross correlation function spectrum 
16384 points were randomly taken and the average am- 
plitude was computed. In order to insure that no data 
points with signal were taken, an iterative procedure re- 
moved all the points which exceeded the average by more 
than 4 times. The average amplitude of noise in the spec- 
trum was considered as an estimate of th e mathematical 
expectation of the noise, £{z) = y/ 7r/2 a. In the case 
when all points of the spectrum are statistically indepen- 
dent^ the cumulative probability for the maximum of n 



aent, the cumulative probability 
points is l|ThomDSon et alJl200lD 



1 



(2) 



where gi is the antenna gain. This factor converts di- 
mensionless fringe amplitude to correlated flux density. 



Then, differentiating this expression, we find the prob- 
ability distribution density of the ratio of the maximal 
fringe amplitude to the averaged noise at points in the 
spectrum with no signal, s = Zmax/ <z>: 

2 n s2 / s2 s n~l 

Pis)^-^se-^{l-e~^) (3) 

It should be emphasized that expression |31 is valid if 
all points of the spectrum are independent. This a priori 
distribution density is an approximation, which is valid 
only partially due to oversampling and deviation of the 
bandpass from a rectangular shape. Using the sample of 
achieved signal to noise ratios, we computed the a poste- 
riori distribution density. We assumed that the a poste- 
riori distribution can be approximated as a function like 
expression 01 with effective parameters <z>qQ and rigff- 
We found parameters of the distribution which minimize 
residuals between the a priori and a posteriori distribu- 
tion densities in the range of signal to noise ratios [3, 6.5] 
by searching through the space of trial values using the 
brute force algorithm. Plots of the a posteriori maximum 
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of the signal to noise ratios and its approximation for the 
case of accumulation periods lengths of 1 s and 0.1 s are 
presented in figures |S1 and El 



2.5 



1.5 



0.5 



—1 1 1 r- 



— 1 1 1 r- 



—1 1 r- 




Fig. 8. — The distribution density of the achieved signal to 
noise ratios (points) and the regression curve (continuous line) 



mulation periods 1.0 s. <2gff>= 1-101, n|,g- : 



"cff-l 



for scans with accu- 
183000. 




Fig. 9. — The distribution density of the achieved signal to 
noise ratios (points) and the regression curve (continuous line) 

p(s) = ^ ^^^^ se IT I 1 ■ 



TT <^eff> 
mulation periods 0.1 s. <2^efT^" 



for scans with accu- 
1.091, rigfj = 989000. 



Using effective parameters of the a priori distribution, 
we can compute the probability that a given observation 
with the signal to noise ratio s belongs to the population 
of observations without signal from a source as 



N{s) = 1 - / p{s)ds = 1 



1 



2 \ Tlott 

1-e--) (4) 



Examining plots of the a posteriori distribution den- 
sity, we conclude that at the range of signal to noise 
ratios [3.5, 6.5] the share of points which belongs to de- 
tected sources is insignificant. Neglecting the population 
of observations with detected sources slightly inflated our 
estimates of the probability that an observation belongs 
to the population of data without signal and makes our 
estimates of the probability of detection a bit more con- 
servative. 



Each source was observed m times, m typically being in 
the range 3-6. At a given series of signal to noise ratios, 
a source has k observations with signal to noise ratio 
greater than sum ■ The probability of false detection pum 
corresponds to it. Then for the case when pum is small 
that we can neglect terms of p^j,„, the probability that 
this may happen due only to the noise in the data, i.e. a 
source was not detected in neither of k observations, is 



i—k 



= Ctl[Nis.O 



(5) 



We used this expression in order to compute the prob- 
ability of false detection for each source. The value 
Piim = 0.05 was used for computations. 

4. RESULTS AND DISCUSSION 

We split the list of 2494 sources into three groups: not 
detected if the probability of false detection is higher 
than 10~^; marginally detected if the probability of false 
detection is in the range of [10~^, 10"^]; and reliably 
detected if this probability is less than lO^'^. 

Careful analysis revealed that 16 detected sources are 
known water masers. The H2O line was within the 
recorded bandwidth. We made a trial fringe search us- 
ing only one frequency channel which included the water 
maser line. The fringe amplitude decreased by factor of 4 
for continuum sources and increased by the same factor 
for maser sources. 

The catalogue of 533 detected sources, which are not 
masers, is listed in tabic |2 Source coordinates are taken 
from the NVSS catalogue (Condon ct al. 1998). Column 
/5/ shows the detection status: D — reliably detected, 
M — marginally detected. The correlated flux density 
in jansky is at three ranges of baseline projections are 
given in columns /6/, /7/, /8/: in the range of [0, 70] 
megawavelengths; in the range of [70, 100] megawave- 
lengths; and in the range of [100, 200] megawavelengths. 
Column /9/ gives the probability of false detection, if it 
exceeds 10~^, or zero, if it is less than 10^^. Column 
/lO/ gives a flag denoting whether the source was previ- 
ously detected in the X/S surveys (KND) or not (NEW). 

Table 131 lists 1945 non-detected sources. Columns /5/, 
/6/ and /7/ provide the minimum, average and maximal 
upper limit of the correlated flux density. These quanti- 
ties are computed on the basis of the calibrated gain of 
the interferometer considering what correlated flux den- 
sity would provide a signal to noise ratio of 6.0. Column 
/8/ gives the number of observations used in the anal- 
ysis, and column /9/ gives a flag denoting whether the 
source was previously detected at the X/S survey (KNO) 
or not (NEW). 

The statistics of the sample of observed sources is given 
in table H VERA detected 60% of the known X/S 
sources and 10% of the new sources. These results are 
somewhat disappointing. With the best weather condi- 
tions VERA is able to detect sources as weak as 100 mJy 
with 2 minute integration times. However, even in win- 
ter, adverse weather conditions significantly deteriorated 
the detection limit. We show in figure [1111 the cumulative 
probability function that a source will have a fringe am- 
plitude signal to ratio limit greater than 6.0 among 14 
observing sessions. This detection limit was computed on 
the basis of the calibrated gain. This probability distri- 
bution characterizes the current averaged VERA single 
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TABLE 2 

List of detected continuum spectrum sources 



jzuuu iitirric 


jjiyou ridiiie 




a (J2000) 


(5 (J2000) 




\P) 


(7\ 
\') 


\P) 


v>) 


{W) 


JOOOl+6051 


2358+605 


00 


01 


07 


09 


+60 


51 


22 


8 


D 




0.16 


0.21 


0.0 


KNO 


J0005+5428 


0002+541 


00 


05 


04 


36 


+54 


28 


24 


9 


D 




0.32 


0.28 


0.0 


KNO 


J0006+5050 


0003+505 


00 


06 


08 


29 


+50 


50 


03 


4 


D 






0.17 


0.0 


NEW 


J0014+6117 


0012+610 


00 


14 


48 


79 


+61 


17 


43 


5 


D 


0.29 


0.25 


0.33 


0.0 


KNO 


J0019+2021 


0017+200 


00 


19 


37 


85 


+20 


21 


45 


6 


D 




0.89 


0.91 


0.0 


KNO 


J0019+7327 


0016+731 


00 


19 


45 


78 


+73 


27 


30 





D 


1.24 


1.35 


1.28 


0.0 


KNO 


J0021-1910 


0018-194 


00 


21 


09 


37 


-19 


10 


21 


3 


M 


0.53 


0.09 




5.D-03 


NEW 


J0027+5958 


0024+597 


00 


27 


03 


28 


+59 


58 


52 


9 


D 


0.35 


0.49 


0.28 


0.0 


KNO 


J0037-2145 


0034-220 


00 


37 


14 


79 


-21 


45 


24 


5 


D 




0.33 


0.12 


5.D-06 


NEW 



Note. — Tabic |^is presented in its entirety in the electronic edition of the Astronomical Journal. 
A portion is shown here for guidance regarding its form and contents. Units of right ascension are 
hours, minutes and seconds; units of declination are degrees, minutes and seconds; units of flux density 
is jansky. 



TABLE 3 
List of not detected sources 



J2000 name 


B1950 name 




a (J2000) 


(5 (J2000) 


(5) 


(6) 


(7) 


(8) 


(9) 


JOOOO+3252 


2358+326 


00 


00:49.72 


+32:52:56.5 


0.11 


0.15 


0.18 


3 


NEW 


JOOOO+5539 


2357+553 


00 


00:20.45 


+55:39:08.6 


0.13 


0.15 


0.16 


3 


NEW 


JOOOl+6443 


2358+644 


00 


01:14.95 


+64:43:01. 1 


0.20 


0.33 


0.88 


9 


NEW 


J0002+5510 


2359+548 


00 


02:00.47 


+55: 10:38.0 


0.15 


0.15 


0.16 


2 


NEW 



Note. — Table HI is presented in its entirety in the electronic edition of the Astronomical 
Journal. A portion is shown here for guidance regarding its form and contents. Units of 
right ascension are hours, minutes and seconds; units of declination are degrees, minutes 
and seconds; unit of flux density is jansky. 



TABLE 4 
Statistics of observed sources 



Category 


Detected 


Not detected 


New, continuum spectrum 


180 


1721 


Known, water masers 


16 




Known X/S, continuum spectrum 


353 


224 


Known X/S/K, amplitude calibrators 


73 


3 


Total 


549 


1945 



beam sensitivity in winter when observations are carried 
out with the integration times of 120 s. 

Since the 2005f_astro X/S catalogue contains 3297 
objects within the dechnation zone > —40°, the zone 
which the VERA can observe, and approximately 60% 
of the sources from this catalogue were detected, we can 
roughly estimate that the VERA can detect 2000 sources 
from that hst. To date, 252 sources were detected at 
K-band with the VLB A in the K/Q survey, 280 other 
known X/S calibrators and 180 new sources were de- 
tected with the VERA in this campaign, and 64 sources 
were detected with VERA at various other experiments. 
In total, the pool of confirmed K band calibrators is 776 
objects. 

5. CONCLUDING REMARKS 

In the 2005/2006 winter campaign ~ 2500 sources were 
observed with the VERA during 148 hours. The target 




Flux density (Jy) 

Fig. 10. — Cumulative probability function for a source to have 
signal to noise ratio above the detection limit 6.0 over 14 observing 
sessions with VERA at 22 GHz. 

sources were either within 6° of the Galactic plane, or 
within 11° of the Galactic center, or within 2.°2 of SiO 
or H2O masers. Among them, 533 continuum spectrum 
and 16 maser sources were detected, including 180 con- 
tinuum spectrum sources not previously observed with 
VLBI. The estimates of the correlated flux densities in 
the ranges [0, 70], [70, 100], and [100, 200] megawave- 
lengths were evaluated. For 1945 non-detected sources 
the upper limits of their correlated flux densities were 
evaluated. 

The distribution of detected sources over the sky is 
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shown in figure ^2 All detected sources are scheduled 
for observations with the VLBA at 22 GHz in June- 
October 2006 for imaging and determination of their po- 
sitions with 1-10 nanoradian accuracies. Tables of de- 
tected and non-detected sources, plots of system tem- 
peratures, fringe plots, and other auxiliary information 
related to this campaign can be found on the Web at 
|tttp : / /lacerta.gsf c .nasa.gov/vlbi/fss 

Analysis of the VERA fringe survey data shows that 
the probability of detecting a 100 mJy source with 120 s 
of integration time is about 10%; the probability of de- 
tecting a 200 mJy source is 60%; and the probabil- 
ity of detecting a 300 mJy source is 80%. Future de- 
velopment a next generation water vapor radiometer 



(Kawaguchi (2006), private communication) promises to 
allow increasing the integration time and, thus, to im- 
prove VERA sensitivity. 



The authors are thankful to D. Gordon for valuable 
comments. The authors made use of the database 

iI] ll997D of the Special Astro- 



CATS (iVerkhodaTir 

physical Observatory. This research has made use of the 
NASA/IPAC Extragalactic Database (NED) which is op- 
erated by the Jet Propulsion Laboratory, California In- 
stitute of Technology, under contract with the National 
Aeronautics and Space Administration. 
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Fig. 11. — The distribution of detected sources over the sky. Disks designate new compact source not observed before with VLBI, circles 
designate known sources observed before at X-band with the VLBA. 



